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We present a study of the decay of an intermediate mini black hole at the first impact of a cosmic ray particle 
with the atmosphere, in the context of D-brane world scenarios with TeV scale gravity and large extra dimensions. 
We model the decay of the black hole using the semiclassical approximation and include the corrections coming 
from energy loss into the bulk. Extensive simulations show that mini black hole events are characterized by 
essentially different multiplicities and wider lateral distributions of the air showers as a function of the energy of 
the incoming primary, as compared to standard events. Implications for their detection and some open issues on 
their possible discovery are also briefly addressed. 



1. INTRODUCTION 

At a cursory look, the inclusive cosmic ray 
spectrum, which has been measured over a wide 
range of energy, seems to be pretty simple in its 
functional shape, characterized by a (fast-falling) 
power-like behaviour over a huge energy range. 
Several mysteries may be hiding behind this ap- 
parent simplicity. Two of these mysteries, which 
have been puzzling theorists and experimentalists 
alike for several decades, and have been exten- 
sively discussed in this Symposium, are the knee 
(10 15 eV) and the ankle (10 19 eV). A long de- 
bate which has spurred from these observations 
has not reached any final conclusions as to the 
origin of the peculiar spectral dependence mea- 
sured in these regions. Physics at the knee may 
involve, according to some interesting proposals, 
a critical behaviour of the QCD pomeron £Q, or, 
alternatively, disoriented chiral condensates and 
production of strangelets or other exotics, while 
the existence or the absence of the GZK cutoff 
P] still needs to be completely assessed. Evading 
the cutoff using string relics [S],0] has also been 
proposed as an interesting possibility. 
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On the other hand, recent formal suggestions 
from string theory and gauge theory [5] lead to 
the exciting possibility of a brane picture of our 
world, with a gravity scale lowered from 10 19 GeV 
down maybe to a fraction of a TeV and with 
the possibility of envisioning exotic new scenarios 
where "black hole resonances" might form even 
at colliders such as the LHC In cosmic ray 
physics, in particular, the energy available at the 
impact of the primaries with the atmosphere can 
be very large and up to several hundreds of TeVs, 
and this points favourably towards a possible test 
of these new scenarios 0. 

We have tested the implications of this hypoth- 
esis for cosmic ray air showers, by studying, in 
particular, the lateral distributions and the mul- 
tiplicities of events mediated by an intermediate 
black hole. Here we provide a brief summary of 
our results, while a detailed analysis can be found 
in®. 

2. MINI BLACK HOLES IN COSMIC 
RAYS 

Various studies of mini black hole production 
at the LHC have been pulished in the last few 
years (see 8 for a rather complete list of ref- 
erences). At the same time analytical and mi- 
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merical studies of greybody factors and the for- 
mation of trapped surfaces in extra dimensional 
models have been able to provide grounds for 
the experimental searches: in the context of the 
brane-world scenario black holes will form copi- 
ously at hadron colliders, in events characterized 
by a large multiparticle multiplicity and a "fire- 
ball" structure, features which would make this 
black hole intermediate state look quite distinct 
from other ordinary resonances. 

Mini black holes form, according to Thome's 
hoop conjecture (see [B] and references therein), 
when the energy in the center of mass of a 
collision is concentrated into a region of radius 
smaller than the corresponding Schwarzschild ra- 
dius. The process may be accompanied by a size- 
able amount of gravitational energy loss, with 
a collision which can be characterized by the 
scattering of two gravitational shock waves of 
Aichelburg-Sexl type at a small impact param- 
eter. The so-formed black hole decays democrati- 
cally mostly on the brane (our world) into all the 
fundamental states of the standard model model 
or of the supersymmetric standard model, if su- 
persymmetry is also open at that energy scale. 

The analysis that we have performed is not 
based on a Monte Carlo modeling of the decay of 
the black hole 9 , rather we write down and com- 
pute the istantaneous parton and lepton emis- 
sions of the decay of the black hole using a multi- 
nomial probability distribution for the hadroniza- 
tion of a given parton, and using suitable frag- 
mentation functions. These are defined via a 
renormalization group evolution up to the mean 
energy of the immediate decay products of the 
black hole evaporation, which is given by the ra- 
tio Mbh/N of the black hole mass Mbh divided 
by the average number N of particles produced in 
such a decay, which in turn is determined in the 
semiclassical approximation by the entropy of the 
black hole. 

3. A MODEL FOR A MINI BLACK 
HOLE DECAY 

As stated above, mini black holes will be pro- 
duced in the collision of any primary cosmic ray 
particle of appropriate energy with a quark or 



gluon of a nucleon in the atmosphere. Such a col- 
lision may lead to the formation of trapped sur- 
faces and of an event horizon, whose size can be 
of the order of 10 -3 fm for a fundamental gravity 
scale M* of the order of a TeV. A spherical s- 
wave forms due to the fragmentation of the black 
hole and the number of states emitted during the 
decay is computed according to a standard for- 
mula for the multiplicities 6 , which depends on 
the number of extra dimensions n. We work in 
D = 4 + n spacetime dimensions. 

We compute the properties of the decay prod- 
ucts, after hadronization, by simulating a se- 
quence of uncorrelated air showers using COR- 
SIKA/SIBYLL |1UI11| and compare the results 
of the multiplicities and the lateral distributions 
of the various components with benchmark sim- 
ulations of standard events. Other hadronization 
models are also available in CORSIKA, for in- 
stance QGSJET [12 . 

4. RESULTS 

Simulations of observables of the type discussed 
in ^31 have been performed at a small cluster. We 
have summarized in 4 figures the most salient re- 
sults of our studies, where we try to distinguish 
between standard (benchmark) events and black 
hole mediated events. We plot in Fig. ^ the to- 
tal particle multiplicity measured at detector level 
(5000 m) assuming an impact at 5500 m (with a 
depth X — 517 g/cm 2 ); a similar simulation at a 
higher altitude of 15000 m (with a depth X = 124 
g/cm 2 ) is shown in Fig. At energies of the or- 
der of 1000 TeV, relevant to the Centauro events, 
the total multiplicity for n — 4 is around 100, for 
a black hole which evaporated 500 meters above 
the detector. A black hole produced closer to the 
detector will obviously lead to fewer particles in 
the shower, an overall picture consistent with the 
Centauro observations. Notice also that the mul- 
tiplicites of black hole events are lower than in the 
standard events roughly by a factor of 10 in the 
"low" energy region, while they are larger also be 
a factor of 10 for energies around the GZK cutoff. 
These differences tend to become less pronounced 
for higher altitudes of the initial impact. 

As we have explained in |SJ, a shortcoming of 
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CORSIKA is that the program is not able to pro- 
vide a neutrino primary, and for this reason we 
have used as benchmarks simulations of proton- 
proton collisions with the same energy and at 
the same altitude of the signal black hole events. 
For this reason, the multiplicities of the bench- 
mark events provide an upper bound for the mul- 
tiplicities of a real event, since hadronic events 
have clearly larger multiplicities compared, for in- 
stance, to neutrino events. As we increase the al- 
titude (Fig. |5J) benchmark events and black hole 
events have similar slopes and intercepts, with 
minor differences. Independently of the altitude, 
black hole events are characterized by wider lat- 
eral distributions (Figs. and compared to 
the benchmark result. This is due to the fireball 
structure of the black hole event, with its origi- 
nal s-wave emission. Our results, of course, rely 
on a semiclassical picture of the particle emission 
by the black hole, given by the Hawking formula, 
and we have corrected for gravity loss in the bulk 
as well. 

We summarize by saying that an analysis of 
lateral distributions and total multiplicities as a 
function of energy show a simple linear behaviour 
(in a log/log scale plot) and measurement of these 
observables and study of their correlation will be 
useful criteria to discern among various models 
for the underlying interactions. 
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Figure 1. Plot of the total particle multiplicity 
as a function of E\ The first interaction is kept 
fixed at 5500 m (517 g/cm 2 ) and the observation 
level is at 5000m (553 g/cm 2 ). We show in the 
same plot the benchmark result (with a proton 
as a primary) and the mini black hole result for 
a varying number of extra-dimensions n. 



trino as primaries in the simulations. As we have 
discussed above, the use of neutrino as a bench- 
mark simulation should only render more re- 
marked the differences between black hole events 
and benchmark events in the characteristics of the 
air showers that have been studied here. 



5. PERSPECTIVES 
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the low gravity scale, which, unfortunately, is not 
a prediction but just a parameter of the theory. 
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Figure 4. As in Fig. [3| but this time the first 
interaction occurs at 15000m (124 g/cm 2 ). 



Figure 2. As in Fig. but this time the first 
interaction occurs at 15000m (124 g/cm 2 ). 
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Figure 3. Plot of the the average distance R of 
the core of the shower of photons as a function of 
E\. The first interaction is kept fixed at 5500 m 
(517 g/cm 2 ). 
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